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Explicit and compact expressions describing the reflection and the transmission of 
a Gaussian beam by anisotropic parallel plates are given. Multiple reflections inside 
the plate are taken into account as well as arbitrary optical axis orientation and 
angle of incidence. 


1. Introduction 

Anisotropic parallel plates are extensively used in ellipsometry.0 To precisely describe such 
experments, it is necessary to take into account internal multiple reflections inside these 
platesH and the Gaussian nature of laser beams.0 However, to the author’s knowledge, no 
general expression of a corresponding Mueller matrix can be found in the literature. 

In Ref. the problem is solved for Gaussian beams but in a particular case: uniaxial 
parallel plate tilted around the optical axis (itself located in the plane of incidence). Al¬ 
though very important results are provided, the formalism introduced by the authors cannot 
be generalized to an arbitrary geometrical configuration, i.e. to a rotating tilted birefringent 
plate with an optical axis not necessarily in the plate interface. Moreover, in this work, the 
calculations were carried out in the direct (a:, y, z) space. This feature has two major impli¬ 
cations: 1) effects related to the beam divergence cannot be studied and more importantly, 
2) the fact that a Jones matrix is only defined under a particular approximation cannot be 
pointed out. 

An adequate formalism to carry out the full calculations is indeed the Fourier optics. 
Theoretical ground foiLthe scalar and vector Fourier optics has been set up some time ago 
in a series of articles. Ba Thanks to this formalism and to the 4x4 matrix method of Ref. |^, 
general and compact expressions describing the transmission and reflection of a Gaussian 
beam by anisotropic parallel plates are provided. This is the topic of the present article. 

This paper is organized as follows. In the first section, useful features of the vector 
Fourier formalism are summarized. This formalism is then used in section ^ to derive a 
general expression for anisotropic parallel plates. Only the paraxial approximation is as¬ 
sumed at this stage. A useful approximation, named here ‘scalar Fourier approximation’, is 
then introduced in section This approximation, implicitly introduced in Ref. |^, provides 
simpler formulas and the possibility to define an extended Mueller matrix for birefringent 
parallel plates. Numerical examples are finally given in section p. 
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2. Vector Fourier optics in the paraxial approximation 

All along this article we shall only be concerned with lossless homogeneous anisotropic media 
and monochromatic Gaussian beam. In this section, we start by considering the propagation 
of a Gaussian beams in isotropic media. The main results obtained in Ref. ^ are recalled 
together with additional information required for the following of the present paper. 

To describe the beam propagation, the direct system axis (x, y, z) is chosen in such way 
that z = k/A: where k is the the Gaussian beam center’s wave vector and k = |k| (see Fig. 
H). The origin z = 0 is taken at the position where the beam size is minimum, i.e. at the 
beam waist position. The position vector will be written r= r_i_ + 2 :z with r_i_ = xx +yy 
and where x , y and z are unit vectors along ox, oy and oz axis respectively. 



Fig. 1. Schematic view of the plane wave propagation in the anisotropic slab. For 
the sake of clarity, some of the inner reflected rays are represented by small arrows. 
The plane of incidence coincides with the plane yz. Symbols a± and b± correspond 
to the four possible propagation directions for the incident angle 9i. The four ele¬ 
mentary transverse walk-offs, Aa± and are indicated. The different vector basis 
used throughout this article are also shown: {x , y , z } for the incidesflgurednt and 
transmitted beams, {x ",y ",z "} for the reflected beam, {xj ,yi ,zi } to perform 
calculations of the birefringence effects. Attached to the two former bases are the 
axis of propagation: z axis for the incident beam, z' axes for the transmitted beam 
and z" axis for the reflected beam. 
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In the paraxial approximation, an electromagnetic scalar field amplitude r_L, z = 0) 
can be expanded according to 


V’( r±,0) = -^11 V'(k_L)exp(-ik_L • r)d^k_L 


( 1 ) 


where the time dependence has been omitted and where k_L = + kyy and |k_L| <C k 

so that kz Rs k(l — |k_Lp/(2/c^)). In Eq. (i. V'(k_L) is the scalar field amplitude in the k_L 
space. To satisfy the paraxial approximation, V'(k±) must be such that the integral has 
appreciable values only for |k_L| <C A:. 

At a given plane z 0, the field is given byQ 


r±,z) = exp(^^(P^ + P^)jjp( r_L,0) 

where, as all along this article, the term exp(ikz) has been omitted and where 

P --i— P - -i— 
ox ay 


Eq. (H) can then be written 


where we dehne for convenience 


r±,2:) = P[ii{k^,z)], 


( 2 ) 


V'(k±,^) = 'i/’(k±)exp 


iz 


k2 


Eor a paraxial Gaussian beam, it is easy to show that 


V'C r_L,0) = 




■ exp 


wl 




where wq is the beam waist, q{z) = z + izR is the complex radius of curvature and zr = kw^ 
is the Rayleigh range. 

Eq. (H) is the plane wave expansion of the vector Gaussian beam. Wave vectors of these 
plane waves are defined by 


kpw = k_L + kzT, « k_L + k{l 



(3) 


in the paraxial approximation. However, Eq. (||) does not take into account the vectorial 
nature of the electromagnetic field. To keep the orthogonality between the electric field, 
magn^ic held and wave vector of each plane waves, one must introduce the six component 
vectorQ 


F{ rx,0) 


(E.{ 

r±,0)\ 

Eyi 

r±,0) 

E,i 

r±,0) 


r±,0) 

By{ 

r±,0) 

\B.i 

r±,0)/ 
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and use the following 


Fourier tra.>.format.„J 


F( ■■x,0) = - 


/ k /c 

V'(k_L) exp(-ik_L • r±)expli-^Gx + i^Gy 




/ “1 \ 
02 
0 

-02 

Ol 

V 0 / 


d^k. 


(4) 


where Gx and Gy are 6x6 matrices which are derived from the Poincare group 
Expressions of these matrices can be found in Ref. It should be noticed that the orthog¬ 
onality between the electric field, magnetic field and wave vector of the plane waves only 
holds in the paraxial approximation, i.e. up to the order C)(|k_i_p/fe^). 

Expression of E( r_L, z) is obtained similarly to Eq. ( 0 ) 




V’(k_L,^;)exp(-zk_L • r±)MQ^Q 


with 
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(5) 


Eocusing on the electric field, Eq. (^) can be reduced to a 3 x 3 matrix equation 


^Ex{ r)\ 

Ey{T^)] = J^[V’(k±,z)M3x3Eo] 

r)/ 

= ^ -ip{'k±, z) exp{-ik+■ r_L)M3x3Eod^k_L 


( 6 ) 


with Eo'^ = (oi, 02 , 0) and 


A , -kyh 


Msxs = 


^F 

— ky kx 


■y^x n\ 

op U \ 

fc2_^2 

1 Ky ^ 

2P“ “4P“ ^ 

=t ^ 0 / 

Notices that the same expression holds for the magnetic field. 


( 7 ) 
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One can verify that, for an electric vector polarized along x (i.e. Eq^ = (1,0,0)), the 
integration of Eq. (^) over leads to the results of Ref. |^. The angular divergence of a 
Gaussian beam therefore generates crossed polarization effects. 


3. Application to anisotropic layer 


We shall now consider a Gaussian beam crossing a single anisotropic parallel plate of thick¬ 
ness d surrounded by a dielectric medium of optical index Ni. The plate is located at 2 : = Zin 
and tilted with respect to k (see Fig. |^). The optical axis has an a priori arbitrary orien¬ 
tation. To compute the transmitted beam, the vector Fourier optics and the 4x4 matrix 
formalism of Ref. |5| will be combined. In doing so, multiple reflections inside the parallel 
plate will be taken into acconnt. 

As indicated in Ref. polarization effects induced by the anisotropic plate can be 
computed by defining an operator acting on every plane waves constituting the Gaussian 
beam. This is further justified since we do only consider here honmgeneous anisotropic 
media. Hence, writing E 4 ( r) for the transmitted beam, one obtainsQ’0 


Et( r ) = exp 



Mt exp 



V'(k±, 0)M3 x3Eo 


( 8 ) 


In Eq. d^, Mt is a 3 x 3 matrix acting on the polarization state of each plane wave and 
z' is an axis parallel to the z axis with 2 ' = 0 at the exit of the plate (see Fig. [^. This axis 
is used to describe the beam propagation after the plate. Inside the plate, the propagation 
of the plane waves is described by the matrix Mt- Let us mention that the position of the 
z' axis along the exit face of the plate can be arbitrarily chosen since it only introduces a 
global phase shift. 

However, in Eq. (^) Mt is determined in the basis {x , y , z }. For consistency with the 
general 4x4 matrix method,Q Mt must be determined in the plane wave polarization basis. 
The polarization vector basis is denoted by {spw, Ppwj kpw} where 


Spw - 


kpw ^ R 


kpw ^ R| 


Ppw — kpw X Spw 


correspond to the TE and TM waves respectively and where n is the unit vector normal to 
the interface. The direction of the wave vector reads kpw = kpw/A: with kpw given by Eq. 


and by convention, {spw,Ppwjkpw} = {x ,y ,z } when k^ = 0. 

The plane of incidence being related to the plane wave vector one gets 


Mt = 9.MtW 


with 
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kpw 
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Mt = 



0 

0 
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(9) 


( 10 ) 


( 11 ) 


and where Ait is the 2x2 matrix describing the transmission of the plane waves in the 
{spw) Ppw} basis. The matrix Ait is obtained by reducing the 4x4 matrix method of Ref. 
^ to a 2 X 2 matrix algebra as described in the following section. 
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A. 2x2 matrix algebra for plane wave propagation inside anisotropic parallel plates 


Because of linearity, the relations between incident, reflected and refracted plane waves 
amplitudes at an interface between two anisotropic media can be written in a 2 x 2 matrix 
form. For the two interfaces I and II of the single layer of Fig. [l|, we thus introduce the 
following 2x2 matrices: 


Rehection 1 - 

a at interface I 

IZla+ 

(12) 

Transmission 1 - 

a at interface I 


(13) 

Rehection a - 

-> 1 at interface I 


(14) 

Transmission a - 

-> 1 at interface I 

Tal.V. 

(15) 

Rehection a — 

> 1 at interface II 

IZal+ 

(16) 

Transmission a — 

> 1 at interface II 

T'al+ 

(17) 


where the + and — subscripts refer to forward and backward beam propagations respec¬ 
tively. TZ and T matrices describe the reflection and transmission of the plane waves. V± 
are the diagonal phase shift matrices 


V± = 


exp(idkaj_ 

0 


n 


0 

exp(idkfe^ 


n 


where and kf,j_ are the four possible wave vectors,l3 inside the anisotropic medium, 
corresponding to the incident wave vector defined in Eq. (|3|). In the above expression, h is 
the unit vector normal to the interface. 

All the matrices of Eqs. (HO) are different in the general case of a biaxial medium. 
They are determined by the electromagnetic held continuity conditions at interfaces I and 
II of Eig. H according to the method of Ref..l3 

Writing and the electric held vectors for an incident plane wave at interface 
I and the transmitted plane wave at interface II, we obtain 


Epw,t — 


%il+'P+^'dla+ + %i+'P+^TZal_V-TZal^Vj^^Ti 


laj_ + • 


E. 




using Eqs. (PHI)- This expression can be simplihed utilising = (1 — A) 


Xit — 7^1+ 


1 — 'P+^TZal_'P-TZal+ 


-1 




(18) 


with Ep^u^^ — Ad^Ep^u j. 

Let us mention that using the same procedure one can also determine the rehected 
electric vector at the interface I of Fig. ||. 


B. Transmitted beam Intensity 


Turning back to the transmission of Gaussian beams, Eq. (^) can be further simplihed 
when one is interested by intensity measurements. For example, if an optical polarization 
component is located upward the anisotropic medium, the out-coming electric held reads 


Eoiit( r) 



exp 



exp 


4 


JMjMaxsEo 
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where 


J = 


J 


0 0 0 


and J is the 2x2 Jones matrix corresponding to this component.u The total intensity 
measured after this component is given by 


lout - 


d r_i_ 


Wn 


(27r)3 
OsxsEo 


.../exp( 


f .{zin + z'){'k\+k'l) 
' exp i - 


k 


0^,3(ki)Eo* 


exp(i(kj_ — k'j_) • r±)d^k±d^k'j_d^ rj_ (19) 


with Osxs = JMtMsxs and where the symbol * stands for the complex conjugate. 


If the matrix J does not depend on the transverse spatial coordinates, the previous 
equation is simplified 


lout - 


27r 


exp 




OsxsEq 


fi^k_L 


( 20 ) 


after integrations over r_L and k(j_ and using 


(52(k_L - k'_i_) = (27r)2 / / exp(i(k_L - k'_i_) • r^)d^ r_L 


where (5^(k_L — k(j_) is the Dirac distribution. 

Eq. (|^) shows explicitly that the total intensity depends only on the waist and not on 
the beam size inside the anisotropic system. This observation has been already made and 
experimentally tested in Ref. 

Up to now only the transmission has been considered. Equivalent expressions can be 
obtained for the reflection by simply replacing by tbe extended Jones matrix describing 
the reflection. 

Eq. (|^) was obtained under the paraxial approximation. This equation can be used as 
it is but, depending on the required accuracy, one can perform further simplifications: the 
plane wave approximation and the Scalar Eourier approximation. The latter is the topic of 
the following section. 


4. Scalar Fourier Approximation 

Although the vector Fourier optics is a useful formalism to describe the Gaussian beam, the 
cross polarization effects ^e indeed very small! (though being observable but essentially 
in extinction experimental). In addition, the birefringence induced by the beam angular 
divergence is also expected to be small, at least for realistic ellipsometry experiments. 
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To a good approximation, one can then assume that all plane waves constituting the 
Gaussian beam have the same wave vector k (i.e. =0). The polarization effects induced 

by a birefringent plate will then be only related to the direction of the Gaussian beam’s 
center. 

This approximation thus amounts to account for the Gaussian nature of the beam only 
in the calculation of the interference pattern of the beam after the plate, as in Ref. After 
the anisotropic plate, the beam is made of a sum of Gaussian beams transversally shifted 
by a distance (the transverse walk-off), induced by successive internal reflections (see Fig. 
1). Using Eqs. (|Tl-0) , the first transmitted beam can be written 


Eit( r') = 


TTWr 


1/2 


-tZR 


exp 


Tn 


alj 


exp 


q{z' + Zin) 

2 q{z'+Zin) 

0 


—ikx 


2 q{z' 


/-ik{y'-Ab 




+ 


eo 


( 21 ) 


with eo^ = ( 01 , 02 ) and where and are the transverse walk-offs. In biaxial media 
there are four different elementary transverse walk-offs 


^ c Aid sin 6*1 cos 01 ^ 

^a± — ^ 


to 


|ka± • n| 


b± 


c Aid sin 01 cos 0 i 
uj Ikfe^ • n| 


where k(jj_ and k;,^ are the wave vectors inside the anisotropic plate corresponding to the 
incident wave vector k, i.e. the center of the Gaussian beam. We should specify that all the 
matrices of Eqs. (12-17) are also determined with respect to the direction of the center of 
the Gaussian beam. 

When the transmission and reflection interface matrices are not diagonal it becomes 
difficult, if not impossible, to write a general formula for the transmitted beam. However, 
using the following property of the Eourier transform 

'*/’( = J^[V’(k_L,^;)] ^ V’( r_L - Ay , 2 ;) = z) exp{iAky)], (22) 

and taking the inverse Fourier transform of Eq. (|2l|), one obtains 


r')] = 


w, 


2 \ 1/2 
0 


1 


27r J 


. {Zin + z')\^\ 
exp( i - - - ) exp 


% 


alj 


exp(i/cyAa+) 0 

0 ex.p{ikyAb^) 
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eo- 


It is therefore possible to sum up all the transmitted beams in the k_i_ space by introducing 
a new set of 2 x 2 matrices: 


W± = 


exp{iAa^ky) 


0 


0 exp{iAh^ky) 

describing the transverse walk-off of the Gaussian beam’s center (see Fig. |l]). 


(23) 




















Following the method introduced in the previous section, one then obtains 


Et(r') = 


with 


2 tt J 


1/2 




exp i 


. {Zin + 2')kj 


Ta_^V^^TZal_'P-Tlal+'P+'^'^la+ + 


k 


-1^ 


exp 


WqKj_ 


Tal+V, Tia, + 


eo 


w\ 




1/2 


2 tt J 


-1 




exp i 


Xzin + z')kl 


exp 


-wlkl 

V 4 






V- 


W+P 

W-V- 


-1 


Mt = % 


al-t 


1 — 'Pj^^TZal_V-TZal+ 


-1 




In the scalar Fourier approximation the usual Stokes vectorsll3 can be defined 

/ l“lP + |02p \ \-^ty\ 


Sin — 


l“ir “ 1021 

Q,\Cl2 “1“ 

\i{aia2 — a\a2)J 


( + \Ety\‘^)d‘^ r_L \ 


Sant = 


ff{\Et^X - \Ety\‘^)d^ 

ff{Et.E;y + Et,Eiy)d^ rx 

H^EtxEiy — E-f.^Ety)d r±J 


(24) 

(25) 

(26) 

(27) 


where Sin and Sout are related to the incident beam and transmitted beam respectively. 
From Eq. {‘M) and following the calculation steps of section ^ one can then determine the 
extended Mueller matrix of the plate Ms such Sout = Ms Sin- 


Ms = 


with 


^ (Pii+Pi2)/2+ 

(pii + pi2)/2 

(p?i+P?2)/2- 

{Ph + pl2)/‘^ 

Pll, 2 l cos 011,21 + 
7*12,22 COS 012,22 

—7*11,21 sin 011,21 — 
\ 7*12,22 sin 012,22 


(p?i+7*ii)/2- 

(7*?2 + Pi2)/2 

(p?i + /*i2)/2- 

(/*ii + p?2)/2 

7*11,21 cos 011,21- 
7*12,22 cos 012,22 

—7*11,21 sin 011,21 + 
7*12,22 sin 012,22 


7*11,12 cos 011,12 + 
P21,22 COS 021,22 

7*11,12 COS 011,12- 
P21,22 COS 021,22 

7*11,22 COS 011,22 + 
Pl2,21 COS 012,21 

—7*11,22 sin 011 , 22 — 
7*12,21 sin 012,21 


7*11,12 sin 011,12+ ^ 

P21,22 sin 021,22 

7*11,12 sin 011,12 — 
P21,22 sin 021,22 

7*11,22 sin 011,22 — 
Pl2,21 sin 012,21 

7*11,22 cos 011,22 — 
7*12,21 cos 012,21 / 


(28) 


Pij,kl 
COS (pij^kl 
sin (pij^^kl 



1 

1 

^^Pijikl 


exp(^ 

'exp 


2 . 
-Wjky 
2 


[mijinh + m*jmki\dky 
j [rriijmli - m*ijmki]dky 


(29) 


(30) 
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and i,j, fc, Z = 1, 2. The matrix elements rriij are defined by 


\m2i m22j 

One feature related to the integral form of Eqs. (|^^) is that pij^ki 7 ^ PijPki- This means 
that one cannot define a Jones matrix (integrated over kj^). 


5. Example: uniaxial parallel plate 


To illustrate our model we shall now consider the very important example of a quarter wave 
Quartz plate (QWP).0 The equivalence between our approach and the results of Ref. is 
formally proved in Appendix I and details concerning the evaluation of Eq. (^) are given 
in Appendix II. 

To study the accuracy of the plane wave and scalar Fourier approximations we shall 
adopt one of the examples of Ref. a Gaussian laser beam of wavelength A = 0.6328 pro. 
and waist tco = 100 /um. For this wavelength, the Quartz ordinary and extraordinary optical 
indices are no = 1.542637 and Ue = 1.551646. 

Since interference effects are sensitive to the plate thickness, we choose to compare 
two realistic components: a first order QWP {d = 87.6010 pm) and a tenth order QWP 
{d = 719.9686 ;um). Finally, two different polar orientations of the optical axis (see Appendix 
II) are chosen 9c = vr/2 (= optical axis in the plane of interface) and 9c = 7r/4. The remaining 
geometric degrees of freedom are the optical axis azimuth (pc and the incidence angle of the 
center of the Gaussian beam 9i. 

The incident beam is assumed to be polarized along ox so that Eq = x (or Sin = 
(1,1,0,0)). As in Fig. ||, the beam crosses a Quartz plate and we shall first consider that 
an intensity measurement is performed after a perfect linear polarizer aligned along ox. 
We shall denote I\\^Gauss and I\\^pw the corresponding intensities computed according to 
Eq. (28) (scalar Fourier approximation) and in the plane wave approximation. We checked 
numerically that results obtained in the scalar Fourier approximation agree with the general 
expression of ^ (to observe noticeable differences one must consider beam waits as small 
as IQ pm which are outside the scope of this article). 

The relative numerical precision of the results presented below has been estimated to be 
of the order of 10“®. This number was determined by checking the energy conservation and 
by looking at the difference between the plane wave and the scalar Fourier approximation 
at normal incidence (they are similar by construction). 

I\\^Gauss is shown as function of pc and 9i in Fig. §(a-b) and ^(a-b) for the two plate 
thickness and the two orientations of the optical axis. 

As expected,! the interference pattern is denser for the tenth order plate (Fig. |(a) and 
(b)) and the intensity is not vr symmetric in pc when 9c < vr/2 and 9i ^ Q (Fig. |2|(b) and 

ia)). 
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(a) Q=tiI2 
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(j)^ / rad 


1*^ order 


Fig. 2. Intensity of a Gaussian beam measured after a first order QWP and a perfect 
linear polarizer: (a) for an optical axis in the plane of interface, (b) for an optical 
axis inclined by vr/d with respect to the plane of interface. The calculations are 
performed using the scalar Fourier approximation and are shown as function of the 
angle of incidence 6 i and the azimuth angle of the optical axis (/>c. 


Fixing (0c = 7r/2,(/)c = 0) and (0c = '^/^,4'c = I\\^Gauss and/|| p^ are plotted as 

function of 0i in Fig. § (a), (b) respectively for the tenth order QWP. In these figures, 
I\\^Gauss is also computed for two beam waists wq = 100/rm and 200/rm. Sizable differences 
appear which increase with the incident angle but decrease when the beam waist increases. 

To quantify the differences the following ratio is computed 

-^11,Gauss ^\\,pw 
I ,Gauss 

and plotted as function of (pc and 0i in Fig. Fig. ^(a-b) and ^(a-b) for the two plate 
thicknesses and the two orientations of the optical axis. One can sees that increases 

with the plate thickness and the angle of incidence. At 0i ~ 0.2 rad and 0c = 7r/2, large 
differences of the order of 10% appear for the tenth QWP. Variations of 5(/||) with (pc are 
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Fig. 3. Same as Fig.2 for a tenth order QWP. 


also sizable. Especially for 6 ^ < vr/4 where interference amplitudes are badly reproduced by 
the plane wave approximation for small values of I||. 


12 

















(b) Q=%IA (j)j.=7r/4 01 / rad 


Fig. 4. Intensity of a Gaussian beam measured after a tenth order QWP and a perfect 
linear polarizer as function of the angle of incidence 9i. Orientations of the QWP 
optical axis are fixed to: (a) 9c = 7r/2 and 4>c = 0, (b) 9c = 7r/4 and cpc = Full 
lines and dashed lines represent the calculations performed using the scalar Fourier 
approximations for wq = 100 m/r and wq = 200 m/i respectively. Dotted lines show 
the calculations performed using the plane wave approximation. 


Another interesting quantity is the degree of circular polarisation when the beam passes 
through a perfect circular polarizer, instead of the linear polarizer of the above example. If 
the polar^r is circular left, using Eq. (28) and the standard Mueller matrices for perfect 
polarizerst3 the beam intensity reads 


II 


Pu + pIi 

2 


Pii,2i sin(()ii^2i 


for Sin = (1) 1^0, 0) and 9c = vr/2. When cjic ~ vr/4, II is small and can be minimized in the 
{9i,4>c) space.cl We present here as function of 9i and (pc for the tenth order QWP in 
Fig. |^(a). 
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Results of Ref. ^ are recovered although rotating thu^late and rotating the polarization, as 
done in this reference, are not strictly equivalent .113’EJ’EJ In Fig. |^(b) the relative difference 
between plane wave and Gaussian beam intensities 

IL,Gauss ^L,pw 
II ^Gauss 

is shown. Large differences corresponding to small values of II are observed. This demon¬ 
strate the necessity to account for the Gaussian nature of the beam in such particular, but 
important, cases. It is to mention that the scalar Fourier approximation and the general 
paraxial calculations are here also in perfect agreement. 




(a) Q=%I2 


s 

u 

u 
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Fig. 5. Relative difference between the quantity of Fig.2 using the plane wave 
approximation and the scalar Fourier approximation. 
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(a) Q=%I2 



Fig. 6. Same a Fig.5 but the for the quantity of Fig.3. 
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10* order QWP, 0 =7t/2 



(b) 


Fig. 7. (a) Inverse of the intensity measured after a tenth order QWP and a perfect 
circular left polarizer. Calculations are performed using the scalar Fourier approx¬ 
imation and are shown as function of the incident angle 6 i and the azimuth angle 
of the optical axis </>c. (b) Relative difference between the quantity shown in the 
top plot calculated using the scalar Fourier approximation and the plane wave ap¬ 
proximation. For these figures, the optical axis is taken in the plane of interface 
(6»c = 7r/2). 


Finally the ratio /9iiP2i//0ii,2i is computed. It increases with 6 i and the plate thickness 
but, even for the tenth order QWP, it does not exceed 0.2 % for any values of (pc and 
9i = 0.2 rad. Assuming piip 2 i ~ 7 * 11 , 21 ) a Jones matrix can thus be defined with an accuracy 
of the order of a few per mill. 

From this study one may conclude that: 

• the scalar Fourier approximation is very accurate with respect to the paraxial approx¬ 
imation when the beam waist is not too small; 

• the plane wave approximation mainly holds for thin Quartz plates, small incident 
angles, large beam waists and optical axes nearly in the interface plane. If these con¬ 
ditions are not fulfilled, an account for the Gaussian nature of the beam is necessary. 

• It is crucial to account for the Gaussian nature of the beam when high performances 
of QWP are foreseen. 
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More generally, accuracies of the various approximations strongly depend on the medium 
birefringence and geometrical configurations. 

6. conclusion 

General expressions describing the transmission and reflection of a Gaussian beam by 
anisotropic parallel plates have been obtained. The vector Fourier optics formalism of Refs. 

is used as general frame work. Multiple reflections inside the anisotropic medium are 
taken into account using a 2 x 2 matrix algebra derived from the general 4x4 matrix 
formalism of RefJ^. The only assumption supplied for these calculations is the paraxial 
approximation. 

To simplify the calculations, a useful approximation introduced in Ref. is then con¬ 
sidered. This approximation consist in taking into account the Gaussian nature of the beam 
only in the interference pattern of the transmitted (or reflected) beam. Birefringence effects 
induced by the beam angular divergence are then neglected. Accuracy of this approxima¬ 
tion was checked for the particular case of uniaxial Quartz parallel plates. For not too small 
beam waists, no noticeable differences were observed compared to the general expression. 

Precision of the plane wave approximation has also been checked using the example 
of Quartz plates. Here noticeable differences were observed. These discrepancies do not 
trivially depend on the geometrical parameters. Roughly, we concluded that they increase 
with the plate thickness and the angle of incidence. In the case of ellipsometry where a 
high purity circular polarization is foreseen, it was shown that an account for the Gaussian 
nature of the beam is necessary. 

As a general remark, accuracies of the various approximations presented in this ar¬ 
ticle decisively depend on the birefringence of the medium, laser wavelength, geometrical 
configuration and type of energy measurements. They must then be checked case by case. 

Interference effects, observed in the variations of the transmitted beam intensity as 
function of the angle of incidence and optical axis azimuth, suggest that a very precise 
calibration of a birefringent plate can be performed without any other optical components. 
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7. Appendix I 

In Ref. 1^, the transmission of a Gaussian beam though a tilted Quartz plate has been 
determined. The configuration is restricted to a tilt axis perpendicular or parallel to the 
optical axis, itself located in the interface. In this appendix we show the equivalence between 
the results of Ref. ^ and the formalism of section |^. 

The elements of the transmission matrix are written 

hotoi exp{-i{kyAo + 5o)) 
l-rl e^p{-2i{kyAo + So)) 

tietel exp{-i{kyAe + Se)) 

l-rj ex.p{-2i{kyAe + Se)) 
rrieo = 0 


rrioo = 

rriee = 
rrioe = 
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where e and o refer to the extraordinary and ordinary wave respectively. The non zero 
extended Mueller matrix elements are given by: 


\Poo\ — 


271 


1 fr 


mn 


rUe 


exp 


‘ exp 


-w^kl 


-wlkl 


dk,, 


dk^ 


(31) 

(32) 


Poo\ee Sin (poo]ee — 2i I 


]exp( 


-Wokl 




dk^ 


= j [^{mee)^{moo) - lft(moo)3‘(mee)] exp^ 

Poo-,ee COS (j)oo-,ee = \ j CXp ^ 


-wlkl 


dk^ 


(33) 


'-w^kl 


J [^(meeWimoo) + 9(moo)9(?7lee)] CXp 


-wlkl 


dki. 


with rriij = ^{rriij) + and where re, Tq, tei, toi, tie, ^lo are the Fresnel coefficients 

for the particular geometric configuration studied here (expressions can be found in Ref. 

ID. 

In Ref. 0 a Jones matrix 


^ ^/2al 
0 


0 




y/ 2 a 2 exp ( —i 


—zarcsm 


Q3 


(34) 


\/“l“2 


is defined and expressions for ai, 02 and 03 are provided (Eqs. (Al 6 a), (A16b) and (Al 6 c) 
in Appendix I of Ref. 

To demonstrate the equivalence between our approach and the one of Ref. § we must 
show that Eqs. (^), (^) and (^) are integral forms of 2ai, 2 a 2 and 203 respectively. In 
order to prove it one just has to expand the integral kernels and then perform the integration 
over ky. Since the scalar Fourier approximation was implicitly assumed in Ref. one finds: 


\po 


=||l+r 


2r:^ 


^ 1 + 2 cos( 4 (/ 93 ;) exp 
-A2 


-2A" 


wf, 


+ 


cos( 2 (/ 33 ;) exp 


2wl 


+ cos( 6 (/Pa;) exp 


-9A^ 

2wl 


I 

I ' ^ 


1 + 2 cos(4(/?a,) exp 


- 2 A 2 \ 

, , /-SA^hl 




+ 2 cos( 8 (^j;)exp 

+ ... 

35 

wi y 

V J\ 


J 

■ ^lo^ol5 

ipx = (5e, <~Px = do and A « 

2 Aof^ 

2 Ae have 


been adopted. Eq. (^) can be further written as a series, and this leads to Eq. (A16a) of 

Ref. |. 

The series expansion of |/3eeP/2 is obviously obtained by replacing and by ry = 
and ty = tigtei respectively in Eq. (^D- This leads to Eq. (A16b) of Ref. 
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In the same way one finally obtains 


Poo\e 


• sin( 


txty 


'oo\ee 


-A2 


- + exp ( jlr^ sm(3<y9a; - ipy) - ry sin(3<y9y - ^Px) j + 


/-A' 


- sin(3[¥53. - (Py]) + j yy sin(5(^y - 3(px) - Vx sin(5v9a; - Scpy) j + 

(? sin(5(^j, - ^ sm(5¥^^ - + 

\ / \' X I y J 


exp 


exp 

-9A2 

2 wl 


■ sm{7ipy - (fx) - - sin(7v9a; - (py) 


+ 


which leads to Eq. (Al6c) of Ref. ||. 

Let us mention that, as stated in section |^, Eq. (34) represents the Jones matrix of a 
quarter wave plate only under the approximation poo,ee ~ PooPee- 


8. Appendix II 

In this appendix, ingredients for the calculation of the double integral of Eq. (|20|) are given. 
Eollowing Ref. 0 we write 

c = sin 9c cos (pc^i + sin 9c sin (pcYi + cos 9cZi 


for the direction of the optical axis inside the Quartz plate, {xj , yi , zj } is the basis 
attached to the Quartz plate (see Eig. |l|). In this basis, the wave vectors of the plane 
wave (see Eq. (^) and Gaussian beam center are given by: 

kpw = kxXi + {ky cos 9 1 + kz sin 0i)yi + {kz cos 9i — ky sin 9i)ii 
k = A:sin0iyi +/ccos6*izi 


with kz ~ k{l — k^/(2A:^))). The plane wave incident and azimuth angles read 


cos 9ipu) = -y cos 9i —A gin ^ tan ( 

/v /C 


'pw 


kz sin 9i + ky cos 9i 
kx 


The ordinary and extraordinary wave vectors corresponding to kpw and c are deter¬ 
mined using the compact expression of Ref. 14. The electric polarization vectors inside the 
plate are determined using the general formula of Ref. 0 . 

Expressions for the interface matrices of Eqs. (0-0) are formally determined thanks 
to the Mapple software package.0 These expressions are much too lengthy to be reproduced 
here. 

The rotation matrix fl (see Eq. (p!o[)) is obtained from 


Spw = As[(fcy cos 01 + ^2 sin 0 i)x —kxCos 9 iy —/cajsin^iz] 

Ppw = Np[{kxkz cos 9i — kxky sin 9i)'k + {kykz cos 9i + k"^ sin 9i + k'^ sin 9i)y 
+(—k_i_^ cos 01 — kykz sin 0 i)z ] 
kpw = Nk[kx^ +kyy +kzz] 


where Ng, Np and are normalization factors. 
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